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We propose that the finite size of the Kondo screening cloud, ^k, can be probed by measuring 
the charge quantization in a one-dimensional system coupled to a small quantum dot. When the 
chemical potential, fi in the system is varied at zero temperature, one should observe charge steps 
whose locations are at values of fi that are controlled by the Kondo effect when the system size L 
is comparable to ^k- We show that, if the standard Kondo model is used, the ratio between the 
widths of the Coulomb blockade valleys with odd or even number of electrons is a universal scaling 
function of ^k/L. If we take into account electron-electron interactions in a single- channel wire, 
this ratio also depends on the parameters of the effective Luttinger model; in addition, the scaling 
is weakly violated by a marginal bulk interaction. For the geometry of a quantum dot embedded 
in a ring, we show that the dependence of the charge steps on a magnetic flux through the ring is 
controlled by the size of the Kondo screening cloud. 

PACS numbers: 72.15.Qm, 73.21. La, 73.23.Hk 



I. INTRODUCTION 



The Kondo effect can be described as the strong renor- 
malization of the exchange coupling between an electron 
gas and a localized spin at low encrgiesii Below some 
characteristic energy scale, known as the Kondo temper- 
ature Tff, a non-trivial many-body state arises in which 
the localized spin forms a singlet with a conduction elec- 
tron. In the recent realizations of the Kondo effect;^ 
a semiconductor quantum dot in the Coulomb blockade 
regime plays the role of a spin S = 1/2 impurity when 
the number of electrons in the dot is odd. The Kondo 
temperature in these systems is a function of the coupling 
between the dot and the leads and can be conveniently 
controlled by gate voltages. As a clear signature of the 
Kondo effect, one observes that the conductance through 
a quantum dot with 5 = 1/2 is a universal scaling func- 
tion G{T/Tk) and reaches the unitary limit G = 2e^/h 
when T <^ Tk for symmetric coupling to the leads.— On 
the other hand, one expects that in a finite system the in- 
frared singularities that give rise to the Kondo effect arc 
cut off not by temperature, but by the level spacing A if 
A T. For a one-dimensional (ID) system with length 
L and characteristic velocity v, the relevant dimension- 
less parameter is /S./Tk ~ ^k/L, where = v/Tk is 
identified with the size of the Kondo screening cloud, the 
mesoscopic sized wave function of the electron that sur- 
rounds and screens the localized spin. This large length 
scale {£^K ~ 0.1 — l^ia for typical values of v and Tk) is 
comparable with the size of currently studied mesoscopic 
devices, which has motivated several proposals that the 
Kondo cloud should manifest itself through finite size ef- 
fects in such systems j^i^ii 

One familiar property of small metallic islands in the 
Coulomb blockade regime is the quantization of charge. 



FIG. 1: Possible experimental setup. Vdw controls the tun- 
neling t' between the small dot (on the left) and the wire and 
Vg varies the chemical potential in the wire. 



Even if the effects of electron-electron repulsion are ne- 
glected or subtracted off, each electron added to the sys- 
tem costs a finite energy because the energy levels are 
discrete. Consequently, the number of electrons changes 
by steps as one varies the chemical potential and sharp 
conductance peaks are observed at the charge degener- 
acy pointsi^ In AlGaAs-GaAs hcterostructures, the level 
spacing of a ID system with size L ~ l/im is of order 
A ^ IQQfieV , large enough to be resolved experimen- 
tally. Clearly, the energy levels of the ID wire should be 
affected by the Kondo interaction with an adjacent dot. 
Therefore, the addition spectrum should exhibit signa- 
tures of the finite size of the Kondo cloud. 

We consider the geometry shown in Fig. [TJ One of 
the ends of a wire of length L with a large number of 
electrons (such that A ^ ej?, where ep is the Fermi en- 
ergy) is weakly coupled to a small quantum dot. The 
dot- wire tunneling is controlled by a gate voltage Vdw 
Both the dot and the wire (which can be thought of as 
a very long and thin dot) are in the Coulomb blockade 
regime. The dot has a very large charging energy and 
a ground state with S = 1/2. We assume that the wire 
is very weakly connected to a reservoir and that trans- 
fer of a single electron between wire and reservoir could 
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be measured by some Coulomb blockade technique. This 
coupling is assumed weak enough so as not to affect the 
energy levels of the wire-dot system. Experiments with 
small dots and large dots have been performed, for ex- 
ample, in Q and considered theoretically in QQ- 

Here we investigate the charge quantization in the 
wire-dot system using the basic Kondo model, ignoring 
charge fluctuations on the dot and assuming the wire con- 
tains a single channel. First we deal with the simplest 
case and ignore electron-electron interactions in the wire. 
In sections [IT] and IIIIl we apply field theory methods to 
calculate the charge steps for the noninteracting case in 
both limits L <^ and L ^k- In section |TVl we 
compute the exact charge steps by solving numerically 
the Bethe Ansatz equations for the Kondo model. This 
approach complements the field theory picture in the in- 
termediate regime, L ^ ^k- The collapse of the numer- 
ical results confirms that the ratio of the widths of the 
charge plateaus with odd or even number of electrons is 
a universal scaling function of £,k/L. In section |Vl we 
include short-range electron-electron interactions in the 
wire using the Luttinger model and discuss how the pre- 
vious results arc modified. In section IVIl we consider 
the same problem for the geometry of a ring coupled to 
an embedded quantum dot, where the flux dependence 
of the charge step locations can also be studied. Section 
IVIIII presents the conclusions. 



II. WEAK COUPLING 

The tunneling between the wire and the small dot in 
the setup of Fig. [T] is well described by the Anderson 
model: 



L-2 



H 



Eh (sVi + '^-^ 
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(1) 



Here the wire is treated as a chain with L ~ 1 sites 
and hopping parameter t. The electron at site j ~ 1 
is coupled to a localized state in the dot with tunneling 
amplitude t' . The parameters Eq and U correspond to 
the energy and the Coulomb repulsion for electrons in 
the dot, respectively. In the Coulomb blockade regime 
t' <^ —Eq ^ Eq + U, the dot is singly occupied and the 
Anderson model is equivalent to the Kondo model^ 



H 



L-2 

E 



Jcl-c,-S, (2) 



where S = d){a/2)d is the spin operator of the electron 
in the dot and J ^ t'"^ / \eq\ > is the antiferromagnetic 
Kondo coupling. In the following we assume that J is 
independent of the chemical potential = eV^, which in 
practice may require that Eq and t' be tuned accordingly. 



For J = 0, the system reduces to an open chain and a 
free spin {S = 1/2). In the continuum limit, we linearize 
the dispersion about the Fermi points and introduce the 
right- and left-moving components of the fermionic field 
"^{x) for electrons in the wire 



* {x) ^ e"=^^V^ij {x) + e-^'^^^^L {x) 



(3) 



The free Hamiltonian for the open chain in the continuum 
limit becomes 

Ho = -ivFj dx (i,'l^dj;'iljj^ + ipldj;ip^^ . (4) 

The open boundary conditions (0) = 5* (L) = im- 
ply that ipL.R are not independent. Instead, -ip^ can be 
regarded as the extension of ipu to the negative- a; axis 



(5) 



We can then work with right movers only and write down 
an effective Kondo model H = Hq + Hx in terms of ijjR 
only (we drop the index R hereafter) 



H 



ivF J dxiP'' d^ip + 2Trv F XqiPHO) T^-ip (0) ■ S, (6) 



where Aq = 2J siii^ kF /t^vf is the dimensionless Kondo 
coupling. For Aq ^ 1, the size of the Kondo screening 
cloud (defined in the thermodynamic limit) is exponen- 
tially large: £,k ^ [vf / D) e^/^° , where D ^ ej? is a 
high-energy cutoff. The boundary condition at the weak 
coupling fixed point [L/^^k 0) reads 

ij{-L) = e^^'ij{L), (7) 

where 5 = kFL mod tt. Using the mode expansion 
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wc obtain the momentum eigenvalues fc„ = (nvr — 5) / L, 
neZ. 

We denote by (/i) the total number of electrons in the 
system, including the one in the quantum dot. We are in- 
terested in the elementary steps of N around some initial 
value Nq = N (/i^)- At T = 0, we calculate (/i) as the 
integer part of N that minimizes the thermodynamic po- 
tential n{N) =E {N) - nN, where E (N) is the ground 
state energy. Defining as the value of /x where A^ 

changes from Aq -I- ^ to Nq + £ + 1, it follows from the 
charge degeneracy condition il {No + £ + 1) = {Nq + £) 
that 



fi^.i ^E{No + £+l)-EiNo + £). 



(9) 



For Ao = 0, we set halfway between the highest occu- 
pied and the lowest unoccupied energy level of the open 
chain. In this case, A^o = odd: the ground state is dou- 
bly degenerate (total spin Stot — 1/2) and consists of one 
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FIG. 2: (Color online) Charge quantization steps for the wire 
coupled to a small quantum dot. The arrows indicate the 
direction in which the single steps move as A (L) grows. 



electron in the dot and pairs of conduction electrons in 
the Fermi sea. The energy levels measured from arc 



(n-l/2)A, 



(10) 



where A = ■nvp / L and n is an integer. The energy of the 
ground state for = A^o + ^ electrons to 0{1/ L) is 



n=l 
{1-\)I2 



ri=l 



+m+i)/2 + 1^*0. (^odd) (11) 

where Eq = E{Nq). Using Eq. dTO]), we get 



E{No + Ao = 0) = + A + s'j + (12) 

where s = for £ even and s = 1/2 for £ odd. It is easy to 
see that, as we raise fi > fiQ, N jumps whenever fi crosses 
an energy level e„. Moreover, there are only double steps 
due to the spin degeneracy of the single-particle states, 
i.e. fJ,2m+i- — M2m+^- ^ result, N (/x) — Nq takes on 
only even values and TV is always odd (see dashed line in 
Fig. HD. 

We obtain the charge steps in the weak coupling limit 
by calculating the correction to the ground state energy 
using perturbation theory in Aq In the limit Ao — > 0, 
the ground state for N odd is still doubly degenerate and 
we write 



\GS{N = odd)) = Yl 
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|0) ® |7) , (13) 



for the values of fc„ defined below Eq. ([5]). The spin state 
of the dot is I7) = Iff) , |-|J-). For N even, there is one 



single electron occupying the state at the Fermi surface 
(with momentum kp). To lowest order in degenerate 
perturbation theory, the ground state for Ao is 



|G5 (TV = even)) = n 414110)® 



(14) 



n<0 



where \s) = Wkp T) ® |-IL) - jfcf i) ® |tr)]/\/2 is the singlet 
state between the spin of the dot and the electron at 
kp- This state has Stot = and { S ■ Sei) = —3/4, 
where Sei = '^If'-fc total spin of the conduction 

electrons. Since the Kondo effect only involves the spin 
sector, the ground state energy must assume the general 
form 



EiNo + £;X) = Eo + j[f 



+ /,(Ao,L)] (15) 



According to Eq. for Aq = 0, /o = and /1/2 = 1. 

For Ao 7^ 0, the singlet formation lowers E {N) for A^ = 
even relatively to A^ = odd. As a result, the Kondo 
interaction splits the double steps and gives rise to small 
plateaus with N even (Fig. To O (Ag), we find 



/(Ao,L) = /i/2(Ao,L)-/o(Ao,i) 

DL 

\ 

Vp 
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(16) 



Note the logarithmic divergence at O (Ag) as L 00, 
characteristic of the Kondo effect. We recognize the ex- 
pansion of the effective coupling A (L) ~ [In (^if /i)]~^ 
in powers of the bare Aq, as expected from scaling ar- 
guments. We have verified that the scaling holds up to 
O (Aq) , i.e., the function / has no dependence on the cut- 
off D but the implicit one in the expansion of A {L) . This 
is remarkable given that E (N) itself is cutoff dependent. 
Based on this, we conjecture that / is a universal scaling 
function of ^k/L. It follows from Eqs. dH), ^ and ^ 
that the charge steps occur at 
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i - li-lfXiL), (17) 



where \_x\ is the floor function or the integer part of x, 
so that 



l£/2\ 



e/2 , i even 
(^-l)/2 , ^odd. 



(18) 



We define the ratio (see Fig. [2]) 

^ EjNo + 3) - 2E{No + 2) + ^(A^o + 1) 
EiNo + 2) - 2E{No + 1) + E(No) 

From Eq. (|16p . we have that in the weak coupling limit 
/(Cx/i » 1) « 1 - 3[ln(eK/i)]-\ (20) 



(19) 
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and the ratio between the width of odd and even steps is The charge steps now appear at 
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L 



> 1 



In 



L 



-ih.ln 



L 



const, 



(21) 

where we included the subleading In hi term in the ex- 
pression of the effective Kondo couphng at scale L. A 
similar scaling function was found in Q for the singlet- 
triplet gap (for fixed N) in the weak coupling limit. 



III. STRONG COUPLING 

When A ^ oo, the spin of the dot forms a singlet with 
one conduction electron and decouples from the other 
electrons in the chain. In this limit the Kondo cloud 
is small (^k/L 0). The strong coupling boundary 
conditions reflect the tt/2 phase shift for the particlc-holc 
symmetric case 



J2S 



(22) 



Note the minus sign relative to Eq. ([7]). This implies 
that the k eigenvalues are shifted with respect to weak 
coupling: fc„ = {ni: + 1: /2 — 5) / L. At the strong coupling 
fixed point we recover the spin degeneracy of the single- 
particle states. The shifted energy levels are 



(23) 



as measured from the original /ig. The ground state en- 
ergy for = iVo + ^ is 



E{No + £) = E[No) + A 



1 



(24) 



where again s = for i even and s 1/2 for ^ odd. 
The ground state consists now of a singlet plus pairs of 
electrons in the wire, thus N (/i) is always even. Fig. [2] 
illustrates how the charge staircase evolves monotonically 
from weak to strong coupling. 

We explore the limit fc]J^^ ^ <C i by working out a 
local Fermi liquid theory^ The idea is that virtual tran- 
sitions of the singlet at a; = induce a local interaction 
in the spin sector of the conduction electrons. The lead- 
ing irrelevant operator that perturbs the strong coupling 
fixed point and respects SU(2) symmetry is 
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(25) 



where the prefactor is fixed such that the impurity sus- 
ceptibility is Ximp = l/(4Tif). This interaction lowers 
the ground state energy when there is an odd number of 
remaining conduction electrons {Sei ~ 1/2, or s = 0), 
thus splitting the charge steps in the strong coupling 
limit. If we introduce the function / as in and 
we find to lowest order in I/Tr- 



-1 + 
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+ 0{{^K/Lf). (26) 
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from which we obtain 
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(27) 



(28) 



At this point we would like to comment about the effect 
of particle-hole symmetry breaking for this wire-dot ge- 
ometry. For the Anderson model of Eq. ^ , particle- hole 
symmetry is absent if the system is away from half-filling. 
We can account for this by adding to the Hamiltonian a 
scattering potential term 



Hp = 2TTVFVij'' {0)ij (0), 



(29) 



where V is of order the bare Kondo coupling Aq ^ l»i 
This term is strictly marginal and can be treat by first- 
order perturbation theory in both weak and strong cou- 
pling limits. Its effect is simply to shift the position of 
the charge steps by 



VA. 



(30) 



Since the shift is independent of the parity of £, the po- 
tential scattering term does not lift the spin degener- 
acy of the charge steps and has no effect on the ratio 
R = Sfi^/fi". 



IV. BETHE ANSATZ RESULTS 

In order to calculate the scaling property of R = 
6 11°/ 5^° for any L/S^k we use the Bethe ansatz (BA) 
solution of the one-channel Kondo problem.— We start 
with a half-filled band of Aq — 1 conduction electrons (A^o 
odd) coupled to a localized impurity spin, corresponding 
to a system size L = (Aq — l)/2 in units where the Bethe 
ansatz cut-off parameter D, related to the bandwidth, is 
set equal to one (see [l3|)- Since the BA solution can 
be obtained for any filling factor, we can add particles 
one by one to the system (A^ = A^o, A^o + 1, A^o + 2, . . .) 
and compute the corresponding energies. This has been 
done by solving numerically the coupled BA equations 
using a standard Newton-Rapson method. In Fig. [3] we 
present results obtained for the ratio R ~ 5^° /Sjj^ with 
Ao== 51, 101, 201, 501, 1001, 2001 and 13 diff'erent val- 
ues of the Kondo coupling in the range 0.06 < J < 0.8. 
The universal scaling curve has been obtained by rescal- 
ing the X-axis L L/^k{J) in order to get the best 
collapse of the data. The entire crossover curve is ob- 
tained, ranging from the strong coupling L 3> to the 
weak coupling regime L <^ S^k- Both strong coupling 
[Eq. [28])] and weak coupling [Eq. ([2T|) ] results are per- 
fectly reproduced. As displayed in the inset (a) of Fig. [31 
BA results agree with Eq. ([^T|) using only one fitting pa- 
rameter: const~ 0.33. The Kondo length scale, shown in 
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FIG. 3: (Color online) Universal ratio R = 5^° /S^"^ as a 
scaling function of S,k/L. Bethe ansatz results obtained for 
various systems sizes {No = 51, 101, 201, 501, 1001, 2001) 
and 13 different values of the Kondo exchange J, indicated by 
different symbols. For each value of J, the system lengths L 
have been rescaled L — > L/^k{J) in order to obtain the best 
collapse of the data using the strong coupling curve Eq. (|28p 
(dashed red line) as a support for the rest of the collapse. 
The weak coupling regime for L <^ £^k, enlarged in the in- 
set (a), is described by the weak coupling expansion Eq. ^21\ 
with constant ~ 0.33 (continuous blue curve). Inset (b): The 
Kondo length scale, extracted from the universal data col- 
lapse of the main panel (black squares), is described by the 
exponential fit Eq. pip with — 0.31 (dashed green line). 



the inset (b) of Fig. [31 displays the expected exponential 
behavior— 



(31) 



with c = 2J/(1 — 3J^/4). Fitting to this expression 
gives ^0 — 0.31 (see inset (b) of Fig. [S]) whieh is in very 
good agreement with the expected value of = I/ti"! 
resultingi. from the impurity susceptibility normalized to 
1/(47^). 



V. LUTTINGER LIQUID EFFECTS 



where V{x) is the screened Coulomb potential and 
Pa{x) = \E'|^(x)\E'„(a;) is the density of electrons with spin 
a. In the low-energy limit, we expand \l/(x) in terms of 
right and left movers as in Eq. ^ and write 

+ [e'''""V^i.(^)^fl.(^) + /i-c.], (34) 

where ipR/La{x) = fpR/Lcr(^)^R/Lcr(^)- Separating the 
processes that involve small momentum transfer (q « 0) 
from the backscattering ones (g « 2kp), we obtain in the 
conventional g-ology notation 



+92\\PB.,aPL,a + — [PR,aPR- 
+92±PR,aPL,-cr} + Hex, 



PL.crPL-a] 

(35) 



where g2\\ = g2± = gi\\ = .94_l = V{0) for a spin- 
independent interaction potential, with V{0) the Fourier 
transform of V{x) at momentum q = 0. The term H^x 
contains the backscattering interaction processes. As we 
will discuss below, this interaction is marginally irrele- 
vant in the sense of the renormalization group. If we 
neglect H^x for a moment, the resulting Hamiltonian is 
exactly solvable by bosonizatioufi^ For open boundary 
conditions, we can work with right movers on a ring with 
size 2L and the boundary conditions of Eq. ([5]). After 
introducing charge and spin densities 



PR,c/s 



PR,] ± PR,l 

V2 



(36) 



one finds that the charge and spin degrees of freedom 
decouple. The effective model for electrons in the wire 
with zero Kondo coupling is the Luttinger model with 
open boundary conditionsii 



Ho + Hi — Hll + Hex , 
where (setting Pq = 0) 



(37) 



q> 
V — c s 



Let us now include the effect of electron-electron inter- 
actions in the wire. We consider the Hamiltonianii 

H = Ho + Hi + Hk. (32) 

The kinetic energy part Hq is given, in the low energy 
limit, by Eq. The short-range interactions are de- 

scribed by 

Hi = IJ2 [ '^^ [ dx' p„ix)V{x-x')p^>ix'), (33) 



where N and S"^ are the zero modes for charge and spin 
excitations, Vc and Vs are the velocities of the collective 
charge and spin modes created by the bosonic operators 
b^ye/s' ^^'-^ is the Luttinger parameter for the charge 
sector. We assume spin SU(2) symmetry and set the 
Luttinger parameter for the spin sector to be Kg = 1. 
From Eq. ([38]) . we have that the energy of the ground 
state for i extra electrons (neglecting Hex^ is 

TT?? 7r7? 

£;(^o+^,Ao = 0) = ii;o + — ^^^ + -^^2, (39) 
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The noninteracting case in Eq. (|T2|) eorresponds to 
Vc = = vp and ~ 1. For repulsive interactions, 
Kc < 1 and Vc > Vg- For weak interactions, we have the 
perturbative resuhs from bosonization 



2^(0) 

TTVp 



-1-1/2 



(40) 



and Vg ~ vp. 

The term H^x on the rhs of Eq. ([57|) is sometimes 
called the marginally irrelevant bulk interaction. This 
is entirely in the spin sector and is analogous to the ex- 
change interaction in the constant interaction model.— 
The operator can be written as 



'^TTgoVs / dxJL{x) ■ Jr{x), 



(41) 



where Jl/b. = ''P^l / / l / r- bare coupling con- 
stant go is proportional to V{2kp) (backward scatter- 
ing process) It is typically very small if the screening 
length is much larger than k'^^ . The renormalized cou- 
pling g{L) obeys the renormalization group (RG) equa- 
tion 



dg 
dl 



-9 



2 



(42) 



where I — ln{DL/vs). Precisely the same interaction 
appears in a spin chain (see, for example, [l3])- The 
solution to 0{g^) is 



9(1) 



90 



i + gol 



(43) 



Therefore, g{L) ~ l/ln(L) vanishes in the limit L — > oo. 
The correction to the energy of the lowest energy state 
due to Hpr ia^ 



SE{No +e,Xa = 0) 



~2Z" 



g{L)s^ + Oig'). (44) 



The Kondo effect has also been studied in Luttinger 
liquids.— One crucial point is that the impurity spin cou- 
ples only to the spin degrees of freedom of the Luttinger 
liquid when the impurity spin is at one end. This is not 
true if the spin couples far from the end. The Kondo 
interaction is as before 



i?K = J*^(£)^*(e)-5, 



(45) 



where e fc^^ <C (e = 1 in the lattice model). In 
general, the spin density operator is 

<i>Hx)^^{x) ~ Jl + Jr+ (e^^'^^^l^^R + h.c}j . (46) 

When we bosonize, the left and right spin densities Jl/r 
can be expressed entirely in terms of the spin boson, but 



the 2kp terms become products of spin and charge op- 
erators. However, for open boundary conditions we can 
use Eq. ([5]) and the entire spin density a,t x = becomes 
proportional to Jr{0). For this reason, spin-charge sep- 
aration is preserved (up to irrelevant operators) for the 
geometry of Fig. 1 with a finite Kondo interaction. The 
bulk Kondo interations only important effects on the spin 
sector, in the low energy limit, are to modify somewhat 
the spin velocity, Vs and to introduce the marginally ir- 
relevant bulk interaction of Eq. (|¥T|) . The finite size 
energies still break up into a sum of spin and charge 
parts. The charge part is exactly as in Eq. In the 

spin part, the dimensionless Kondo coupling is defined as 
Ao = 2Jsm^{kpe)/TTVs, with Vs replacing vp in contrast 
with Eq. The Kondo coupling is still marginally 

relevant for A% = 1. We can then write 



^(A^o + Ao) = 



AR'L 



^,UX,g,L). (47) 



The scaling functions ,fs{x), s = 0, 1/2, are the same ones 
that we calculated ignoring Luttinger liquid interactions, 
apart from the contribution from the marginally irrele- 
vant bulk interaction. To lowest order in A and g, we 
have 



f{\,g,L) = h/2i\,g,L)- foiX,g,L) 
= l-|g(L)-3A(i). 



(48) 



Note that both the bulk marginal operator and the 
Kondo interaction reduce the energy when N is even. 

The RG equation for the Kondo coupling A is modified 
by the marginal bulk operator— 



d\ ^, 



(49) 



Solving this equation in the presence of the second term, 
one finds 



X{L) 



21n(L/Li) 



h-i'iU/Li) - \n\L/Li) + 2ln{U/Liy 



(50) 



where is defined such that A(L = ^k) = 1; and Li 
is a characteristic length scale for the bulk marginal in- 
teraction defined by g{L) « l/ln(L/Li). If the bare go 
is small, then Li ~ {vs/ D)e~'^/s° < Vs/D. In the limit 
Li — > we recover the scaling function A = \{S,k / L) ~ 
[ln(^if/L)]~^. In general, we expect A = X{S,k / L, g{L)). 
Eq. (|49p also implies an unusual dependence of S,k on 
the bare Kondo couplin g^^i^^ 



j\ exp 




where c is a positive constant of 0(1) 
VajD we recover the usual result 



\DLi 

(51) 

In the limit Li <C 
{vs/D)e-^^^°. In 



the limit A — *■ with g held fixed this becomes oc 



^const/ V Ao 



7 



In the charge staircase (see again Fig. [2]), the width of 
the plateaus is normahzed by the charge addition energy 
Ac ~ TTVc/KcL. The charge steps for ^ L arc given 
by 



Mo 



Ac 
where 



l~uf[^,giL) 



VsKc/Vc < 1 



(52) 



(53) 



and / is given by Eq. (|l8|) . As the charge and spin 
addition energies are different for u ^ Iji^ there arc no 
double steps in the interacting case even for A ^ 0. The 
ratio between odd and even steps including interaction is 



R 



l + uf{U/L,g{L)) 
l-uf{U/L,g{L)y 



(54) 



Thus we see that the main effect of the screened bulk 
Coulomb interactions is to suppress R somewhat due to 
the reduction of the parameter u from 1. As m — > (for 
very strong Coulomb interactions) the even-odd effect 
disappears entirely and we see simple Coulomb block- 
ade type steps. For cleaved edge overgrowth quantum 
wires, u « q 5, i9i20,2i carbon nanotubes, one finds, 
typically, u k, Q.l^ (For further discussion of the exper- 
imental possibilities see Sec. VII.) For ^ L, 



R 



l + u- 3uX{L) - lug{L) 
1-U + 3uX{L) + \ug{L) ' 



(55) 



Now consider the strong coupling limit, ^ L. The 
Fermi liquid interaction is the same as in Eq. (j25p . but 
with given by Eq. ([?T|) . In contrast with Eq. pS)) . 
the marginal operator now lowers the ground state energy 
when N is odd, because in the strong coupling limit the 
number of free electrons in the wire in N — 2. Thus we 
have in this case 



fi^K/L,giL)) 



-1 



-ML) 



2L ' 



(56) 



The charge steps for £,k ^ L are given by Eq. ([5^ with 
/ taken from Eq. ([56|) . The ratio becomes 



R = 



1-u 



ug{L) 



' 2L 



1 + u- §ug{L) - u^^ 



(57) 



If we set g ~ 0, we can express the ratio R in terms of 
the ratio calculated for the noninteracting case 



R = 



1 - u + (1 + u)R 
1 + M + (1 - u)R' 



(58) 



where R = R{S,k /L) is the scaling function shown in Fig. 
[31 For u ^ 1, the function R is finite at both weak and 
strong coupling fixed points. In the weak coupling limit. 
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FIG. 4: (Color online) Effect of the marginally irrelevant 
bulk operator on the scaling function f{^K/L,g{L)) calcu- 
lated from the Bethe ansatz. The symbols correspond to the 
noninteracting case, g{L) — 0, as in Fig. O The black circles 
are obtained with the Bethe ansatz solution of the spin chain 
Kondo model ^ ior L = 1000 (p(lOOO) ~ 0.115). 



R ^ oo and i? ^ (1 + u)/(l — u). In the strong coupling 
limit, i? ^ and ^ ^ (1 - u)/(l + u). 

A nonzero go leads to a weak violation of scaling be- 
cause / is not only a function of £,k/L but also of g{L). 
The effective g{L) is typically small in quantum wires, 
except at very low electron densities, and could in prac- 
tice be treated as a fitting parameter. We can write down 
a general expansion for / in powers of g: 

f{iK/L,g) = h{iK/L)+gh{^K/L) + 0{g^). (59) 



It follows from Eqs. (gS]) and ^ that /i(^k/L > 1) 
-3/2 + 0{\{L)) and h{U/L « 1) = 3/2 + OiU/L). 



FIG. 5: Schematic picture for an open Heisenberg chain cou- 
pled to a spin impurity (arrow) at the left boundary via a 
weak antiferromagnetic exchange (dashed bond). 

In fact, the function f{^K/L,g) can be calculated in 
the framework of a Heisenberg spin chain model with a 
weak coupling at the end of the chai n^^i^^ 



L-l 



H = Ji ^ S i-S, 



i+r 



'Hi. 



Jx^imp-Si. (60) 



i=l 



This model is depicted in Fig. [5l Its low energy limit is 
the same as the spin sector of the Luttinger liquid J^i^ii^ 
In this case the bulk marginal coupling has a bare value, 
go, of 0(1). The Kondo coupling is proportional to J'x'^ 
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Since the charge sector separates from the spin sector 
anyway, in the Luttinger hquid Kondo model, it follows 
that the needed function f{£^K/L,g{L)) can be deter- 
mined from the spin chain model. 

The fmrction f (£^k I L , g{L)) can be evaluated via the 
energy difference between the ground states for even and 
odd chain length, of total spin and 1/2: 



AL 



(61) 



One can extract this quantity exactly using the 
BA solutionis of the model ([60|) by computing 
i [E{L + 2) + E{L) - 2E{L + 1)] with L even for various 
Kondo coupling strengths J^. This has been achieved for 
L ~ 1000 as shown in Fig. IH where one can see the full 
scaling function with the two limiting cases 



/ 



1 

-1 



K 



if Lii 

if L/U 







(62) 



In order to obtain the scaling function / of the spin chain 
Kondo model in the entire regime of L/S^k shown in Fig. [4] 
with L = 1000, we had to convert — > using the 
unusual exponential square root behavior for the Kondo 
length scale (jSip . In order to make things quantitative, 
we used the following convcrsionii for £,k{J'k) 



£.K = Co exp 



(63) 



witbi^iiS Co 



It is reassuring to note that, even 



for this large value of bare coupling g^, the effects of the 
bulk marginal operator on the cross-over function / are 
fairly small for L — 1000. We might expect these effects 
to be even smaller for the Luttinger model with realistic 
parameters chosen to describe the systems discussed in 
Sec. VII. 



VI. FLUX DEPENDENCE IN QUANTUM 
RINGS 



where e ~ kp^ <C L. In the symmetric case it is conve- 
nient to label the eigenstate of the open chain by their 
symmetry under the parity transformation x ^ L — x. 
Expanding '^[x) as in Eq. ([3]), we have 

+ e-'^^^VL(£)].(65) 



*(e) + *(L-e) =e" 

I AkpiL — e) 



"''"'[V'/?(£)+e 
[i^R{L - e) 

This motivates defining the even and odd fields 



(66) 



This way, for zero flux Hk involves only the even channel 
since 



^'(e)-h*(L-e) 

71 



V'e(e) + e*''"(^-^Ve(i-e). 



(67) 

For J = 0, the system is equivalent to an open chain 
and a free spin. In terms of even and odd fields, the free 
Hamiltonian of Eq. (|4|) is simply^ 

i?o = -^VF I dx {^Id^i^, + i^ld.,^o) ■ (68) 

^0 

The open boundary conditions in the weak coupling limit, 
*(0) = ^'(X) = 0, imply 

V'e/o(i) = Te-'"-^Ve/o(0). (69) 

Because of the relative minus sign in Eq. the even 

and odd channels are nondegenerate. The energy levels 
alternate between the even and channels. We assume 
that the lowest energy single-particle state in the band 
is even under parity. In this case, for A^o = Ap + 1, p 
integer, the lowest unoccupied state belongs to the even 
channel. The energy levels relative to fiQ for J = are 
el = (2n-hl/2)Aand e° = (271-1/2) A, with A = ttvf/L 
and n integer. 

We can account for the magnetic flux through the ring 
by introducing a phase in the hopping between the quan- 
tum dot and the two leads This does not affect the 
Hamiltonian of the open chain in ([55]) . but the Kondo 
interaction is modified to 



The effect of the Kondo screening cloud on the charge 
steps of a ID system is not unique to the wire geome- 
try. One interesting alternative is to suppose that the 
quantum dot is embedded in a ring with circumference 
L (inset of Fig. [S]). This geometry offers the possibility 
of looking at the dependence of the charge steps on the 
magnetic flux $ threading the ring, which is related to 
the persistent current for the embedded quantum doti^i^ 
In this section we focus again on the nonintcracting case. 
For simplicity, we assume that the coupling between the 
dot and the leads is parity symmetric. We can write the 
Hamiltonian (for zero flux) as H = Hq + Hk , where Hq 
is the same as in Eq. ^ and the Kondo interaction is 



Hk = J[*^(e) + - e)]-[^{e) + *(i - e)] • S, (64) 



o«a!lT;t 



(70) 

where a — 27r$/$o and $o ~ hc/e is the flux quantum. 
The Kondo interaction in terms of even and odd fields 
reads 



cos -iPl (0) 



*sin-i/;J (0) 



a 
"2 



cos-7Ae(0) + isin-V„(0) -S, (71) 



where Ao = 4 J sin (kpe) /ttvp- If particle-hole symmetry 
is broken we have to add the potential scattering term 



H 



p 



2hvfV 



cos-^A^ (0)-i sin -7/^^(0) 



cos -V-e (0) 



(X 

i sin -V'o (0) 



(72) 
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FIG. 6: (Color online) Flux dependence of the charge steps 
in the ring with an embedded quantum dot (inset), assuming 
no potential scattering {V = 0). Dashed lines: weak coupling 
regime {^k S> L); solid lines: strong coupling regime [S^k ^ 
L). 



where V ^ 0{Xq) is the dimcnsionlcss coupling constant. 
This operator is strictly marginal in the noninteracting 



1 and treat Hp using first 



case. We will assume V ^ 
order perturbation theory. 

For A = 0, the energy levels are flux-independent, due 
to the open boundary conditions. In the weak coupling 
limit A(L) ^ 1, perturbation theory in the Kondo inter- 
action of Eq. ([7T|) for general a yields a different splitting 
for even and odd channels. As a result, the spin part of 
the ground state energy depends on £ mod 4. To 0(Ao), 
we find 



E{No + e) = 



Eo + fi*a( + A l^- + . 

~3s^ A {1- cos ae)X{L) 

-3s2(ln2) Acosa^ (1 - cosa^) X^{L) 

+ ^{l + cos'a) A2(L) + 0(a3), (73) 
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where ag = a + {£ + 1)tt/2; s for £ even and s ~ 
1/2 for i odd. In Eq. ^ it is implicit that a Aj^lnL 
term is generated by perturbation theory with the correct 
coefficient to recover the expansion of A(L). 

The charge step locations can be calculated using Eq. 
([9]) . The last term in Eq. ([73]) is independent of i and 
drops out when taking the difference E{No + £ + 1) — 
E{No + £). To first order, the potential scattering term in 
Eq. ((72)) simply increases the energy by yA(l -l-cos a) for 
each electron added in the even channel and by T^A(1 — 
cos a) for each electron added in the odd channel. The 
charge steps in the weak coupling limit are located at 



£ 


1 


1 -1- cos + n 


£ 










)] 


2 


+ 2 + 




_2_ 





V-'-i-lY\iL) + 0{X') 



(74) 



Note that, although potential scattering gives some weak 
flux dependence to the energy levels, it does not lift their 
spin degeneracy. Only the term due to the Kondo inter- 
action alternates between £ even and £ odd and therefore 
lifts the spin degeneracy. Hence the main effect of the 
Kondo interaction in the weak coupling limit is to split 
the charge steps that are degenerate for A = 0, namely 
1^1+1/2 and /i£+3/2 with £ even. We define 



n+i 



(75) 



where n is an integer. From Eq. (|74p we obtain 



2A 



[l + (-l)"cosa]A(L) + Ci(A2). (76) 



For symmetric coupling to the leads only the even chan- 
nel couples to the impurity when a = 2rmT, m integer. 
As a result, the steps in the odd channel (corresponding 
to odd values of n in Eq. ([76)) ) do not split. The opposite 
happens when a = {2m + l)7r (see Fig. [5]). 

In the strong coupling limit, the tt/2 phase shift modi- 
fies the boundary conditions in such a way that the trans- 
mission across the dot becomes perfect For zero 
flux, the phase shift in the even channel implies 



(77) 



With these new, periodic boundary conditions, the right- 
and left-moving components are decoupled 



-ikp L 



V'i?/L(0). 



(78) 



The system is then equivalent to an ideal ID ring, whose 
eigenstates can be labeled as right or left movers. The 
degeneracy of right and left-moving energy levels at flux 
TOTT implies charge steps of magnitude 4 in the strong 
coupling limit, as can be seen in Fig. [5| As we did in the 
weak coupling limit, we can introduce the flux in the hop- 
ping from X = L to X = Q. Since the ring is now closed, 
the unperturbed energy levels are fiux dependent ^21 



E^/\a) = {2n-ll2±a/ir)A, 



(79) 



and exhibit zigzag lines with level crossings at a = tott, 
m integer. The energy levels in Eq. (|79p are measured 
from /Zq defined in the weak coupling limit. The shift of 
— 1/2 is due to the phase shift of 7r/2 for the states in 
the even channel at a = and can be understood in the 
following way. Consider the lattice model with N— 1 sites 
and a symmetric Kondo coupling between the electrons 
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at the ends of the chain, at j = 1 and j = N — 1, and to 
the impurity spin at .7=0, 



L-2 



J I J 



,+c^_, -S. (80) 



For J — and at half-fiUing, there are A^ — 1 free conduc- 
tion electrons in the chain. The hopping Hamiltonian is 
invariant under the particle-hole transformation 



(81) 



For J 7^ 0, particle- hole symmetry is broken by the 
Kondo interaction if N is odd and preserved if N is even, 
since under particle-hole transformation 



Ci -f c 



JV-l 



(82) 



As we chose to fix /Iq so that A'q = 4p + 1 is odd, the 
charge steps in the strong coupling limit are not symmet- 
ric about (jLq. However, the charge steps must be symmet- 
ric about ^Q± A/2, which correspond to values of chem- 
ical potential with an even number of electrons in the 
weak coupling limit. Notice that, if we start from A^o odd 
in the weak coupling limit, there is a value of A(L) ^ 0{1) 
(for L ~ ^k) But which the charge step /ii crosses fiQ for 
a ~ 0. This means that, for fixed fi ~ iiq and a ~ 0, 
the ground state at the strong coupling fixed point has 
an even number of electrons {N ^ Nq + 1 — Ap + 2, i.e. 
A^o — 1 = 4p free electrons in the ring plus two in the 
singlet). 

Using the strong coupling boundary conditions, the lo- 
cal Fermi liquid interaction analogous to Eq. ([^S]) can be 
written 



H 



FL 



II Up 



(83) 



The main effect of this interaction is to lift the degen- 
eracy between right and left movers at a = rmr and to 
split two out of the four values of We calculate 

the correction to the ground state energy using degen- 
erate perturbation theory to 0{1/Tk) for V = 0. The 
calculation is similar to the weak coupling limit of the 
side-coupled quantum dot^ and is presented in detail in 
Appendix A. The charge steps (as labeled in the weak 
coupling limit) for a w rmr, m even, are given by (see 



Fig. El). 

M4ri+l±i ^ f^O 



1 TT^R- 




X - '"+2^V^ + IT 




For a ~ mTT, m odd, the charge steps are located at 




A = 2"+2±^^lL 




Note that two charge steps remain degenerate at each 
point of level crossing a = tott. This is because for a = 
imr either the single-particle states in the odd channel 
(for m even) or the states in the even channel (for m 
odd) decouple from the spin of the dot. As a result, the 
Fermi liquid interaction reduces to 



2tt^ vl 



(88) 



where the index e applies to m even and the index o 
to m odd. The two degenerate charge step correspond to 
adding electrons to the state whose energy is not modified 
by the Fermi liquid interaction. 

If the potential scattering is nonzero, the degener- 
acy between the even and odd channel at a = mi: is 
lifted. For ^k/L ~ V <^ 1, both the Fermi hquid in- 
teraction and the potential scattering can be treated us- 
ing degenerate perturbation theory near a w ttitt. If 
^k/L <^ V ^ 0(1), the Fermi hquid interaction can be 
added on top of the unperturbed energy levels of an ideal 
ring with a delta function potential at a; = 0. In any case, 
the result is that, due to avoided level crossings, potential 
scattering suppresses the flux dependence of the charge 
steps in the limit S,k / L ^ 1. However, a feature that 
survives when 7^ is that the splitting of the double 
steps in the even (odd) channel for a ~ rmr with m even 
(to odd) is proportional to £,k/L. 

This significant difference in the flux dependence of the 
charge steps, between the weak and strong Kondo cou- 
pling limits, is consistent with previous results for the 
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persistent current. The persistent current can be calcu- 
lated by taking the derivative of the ground state energy 
with respect to the magnetic flux at fixed electron num- 
ber: 



3 = -e{dE/da)\N, 



(89) 



where e is the electron charge. Using Eq. ([73]), we recover 
the results derived in Ref. [5[ for the weak coupling limit. 
Here we shall be concerned with the strong coupling 
limit. At the strong coupling fixed point, S,k/L 0, 
the persistent current has a saw-tooth shape as expected 
for an ideal ring (Fig. [7]). For small but finite ^k/L, we 
expect that the Fermi liquid interaction rounds off the 
sharp features near the degeneracy points. 

Consider first the case of N odd. For both mod 4 = 
1 and N mod 4 = 3 the ground state is quasi-degenerate 
for a « mTT, m integer. From first order perturbation 
theory in the Fermi liquid interaction, we find that the 
ground state energy for a w tott and S^k/L <C 1 is (sec 
Appendix A) 



A 




const, 



AL 



(90) 



where [9] denotes the value of 6 reduced to the interval 
— TT < [9] < TT by subtracting an integer multiple of 27r, 
i.e., [61] = 61 if 16*1 < TT, [61] = 61 - 27r if TT < 61 < 37r, etc. 
The persistent current for N odd and a ~ ttitt is 



Jo(a) = — 



2 (a — 771 tt) 



(a — rmT)/TT 



4L 



(91) 

In the above case the persistent current is paramagnetic 
and has tt periodicityi^ 

Now consider A^ even. For A^ mod 4 = 0, the ground 
state is quasi-degenerate for a « rmr, m even. The 
ground state energy near these points is 



E{a) 




const.(92) 



For A^ mod 4 = 2, the degeneracy points appear at a 
7777r, with 777 odd. The ground state energy is 



A 




K 



2 1 a — 7777r I 



4L 



K 



8L 



const. (93) 




FIG. 7: (Color online) Persistent current for the embedded 
quantum dot with A'' electrons, A'' mod 4 = 0, in the strong 
coupling limit. The dashed line represents the saw-tooth 
shape expected for periodic boundary conditions at the strong 
coupling fixed point. The solid line is the result of Eq. (|94p 
for L/C A- = 20. 



Eq. (|93p can be obtained by a translation a ^ a + tt 
of Eq. ([92]) . In both cases, the persistent current for A 
even near a w rmr [m even for A" mod 4 = 2 and m odd 
for A^ mod 4 = 2) reads 



3e{a) = — — 



2(a — r777r) 



2[a — 77i7r)/7r 

\ TT I ^ \ ?,L ) 



The persistent current for A^ even has 27r periodicity. It 
is paramagnetic for A^ mod 4 = and diamagnetic for 
A^ mod 4 = 2. The persistent current for the embed- 
ded quantum dot with N = Ap in the limit <C i is 
illustrated in Fig. \7\ 

The effect of electron-electron interactions on the 
transport through an embedded quantum dot was dis- 
cussed in j5ll28|. For ATc ~ 1, the system still flows 
to the one-channel Kondo fixed point, with one elec- 
tron in the even channel screening the impurity spin. 
However, one important difference from the noninteract- 
ing case is that, for Kc < 1, the backscattering term 
(-- iIj'1{0)iPr{0)) of Hp in Eq. ((721) is a relevant pertur- 
bation in the sense of the renormalization group. As in 
the Kane-Fisher problemf2£ the effective potential scat- 
tering grows with L as 14// fx l/i(i-^^<=)/2 for Veff <^ D 
and as K// oc yi(i/^fc-i)/2 y^^^ ^ r^^ie large 
Vef f suppresses the flux dependence of the charge steps 
in the strong coupling limit of the Kondo effect and the 
persistent current decreases with L faster than 
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VII. EXPERIMENTAL POSSIBILITIES 

With current technology, it is not clear whether one 
can experimentally realize a ballistic single-channel ring 
of a suitable length to directly test the theoretical pre- 
dictions for this geometry. However, prospects are more 
encouraging for the linear geometry, where a straight wire 
segment is coupled to a localized spin at one end. Candi- 
date systems may include carbon nanotubes and wires 
made from semiconductor nanocrystals, as well one- 
dimensional conducting channels fabricated from two- 
dimensional semiconductor structures i^i^i^i^^ Perhaps 
the most promising candidates are one-dimensional wires 
constructed, using cleaved edge overgrowth, at the edge 
of a GaAs single-well or double-well structure.— i^i^iiSS 
Electron mean-free paths of the order of 10 microme- 
ters have been obtained in such structures, at electron 
densities in the range of 50 electrons per micron. Local 
electron densities in the wire can be controlled by means 
of top gates, and one can thereby produce a series of con- 
ducting segments, of various lengths, separated by barri- 
ers of depleted regions. Electrical contact to conducting 
wire segments can be made through a two-dimensional 
electron gas, which is present within the GaAs well in 
regions where it is not depleted by a top gate. In dou- 
ble well structures, one can also establish momentum- 
controlled tunneling contacts between parallel wires in 
the two wells 

In order to realize the geometry of interest here, one 
might create a short conducting segment with a small odd 
number of electrons (a Kondo dot), separated by a de- 
pleted barrier region from a much longer conducting seg- 
ment, which would constitute the one-dimensional wire. 
The strength of the Kondo coupling could be controlled 
by varying the voltage on the top gate above the barrier 
region. Electron energy levels could be studied by tun- 
neling electrons into the conducting segment from a lead 
at the opposite end from the Kondo dot, from a weakly 
coupled two-dimensional gas at the side, and/or by tun- 
neling from a second parallel wire. Tunneling measure- 
ments can give information about excited states of the 
wire system as well as the ground state. Alternatively, 
the charge state of a conducting segment can be moni- 
tored through its Coulomb interaction with a nearby sin- 
gle electron transistor or quantum point contact device 
whose electrical conductance is sensitive to small changes 
in the electrostatic potentiali^ 

Single-wall carbon nanotubes are another promising 
candidate for realizing effects analogous to those dis- 
cussed in this paper. Electron mean-free-paths in nan- 
otubes can be very long. Gated quantum dots have been 
fabricated in carbon nanotubes, '^^i^'^i'^^ and Kondo-type 
effects have been observedi^ Analysis of effects in a car- 
bon nanotube is made more complicated, however, by 
the existence of two orbital conducting channels, which 
arc degenerate in an ideal nanotube, but may be split by 
distortion or imperfections in actual nanotubesi^i^ 



VIII. CONCLUSION 

We have shown that the width of the charge steps in a 
mesoscopic device depends on the finite size of the Kondo 
cloud. In general, the signature of this effect is the broad- 
ening of the Coulomb blockade valleys with total N even 
as we increase Tk (decrease £,k)- For the case of a single- 
channel noninteracting wire, the crossover between weak 
{L <C £,k) and strong {L ^ £_k) couphng is described by 
a universal scaling function R{^k/L), defined as the ra- 
tio between the width of the odd steps and the width of 
the even steps. Using perturbation theory in the Kondo 
interaction in the weak coupling limit and perturbation 
theory in the effective Fermi liquid interaction in the 
strong coupling limit, we derived the asymptotic behav- 
ior of R{S^k/L) for L and L £^k- These formulas 
are in agreement with the exact numerical results ob- 
tained using the Bethe ansatz solution. We generalized 
these results for the case of a quantum dot coupled to a 
Luttinger liquid. If we neglect the effect of the marginal 
bulk interaction, the ratio R between even and odd steps 
is still a scaling function of S^k/L, but the asymptotic 
values of R for L <^ and L ^ are determined by 
the parameter u = VsKc/vc- More generally, R is also 
a function of the effective coupling constant g{L) associ- 
ated with the marginal bulk interaction. Finally, we have 
shown that, for the geometry of a quantum dot embed- 
ded in a mesoscopic ring, the charge steps are weakly fiux 
dependent for L <^ and strongly flux dependent for 
L S^K- The latter behavior is reflected in the persis- 
tent current near the strong coupling fixed point, which 
we calculated using a local Fermi liquid theory. 
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APPENDIX A: CHARGE STEPS FOR 
EMBEDDED QUANTUM DOT IN THE STRONG 
COUPLING LIMIT 

In this appendix we show how the Fermi liquid interac- 
tion of Eq. ((83l) lifts the fourfold degeneracy of the charge 
steps of the ideal ring near the points of level crossing. 
We focus on the calculation of groimd state energy E{N) 
for fixed A'^ and a ~ using degenerate perturbation the- 
ory. We generalize the results for a « tott, m integer, at 
the end of this section. 

We start by separating the mode corresponding to the 
highest partially occupied energy level from the remain- 
ing ones that are either completely occupied or com- 
pletely empty. We write A^ = A^o -I- = 4(p -I- n) -)- 2 -|- 
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with n integer and £' = {N — 2) mod 4. Then the 
quasi-degenerate levels at the Fermi level are en^^(a) = 
2A(n+3/4±a/27r) > for > iVo + 1 (see Fig. El- (Here 
we measure the energy levels from /Iq and set ^j,q — for 
a shorthand notation.) We rewrite the Fermi liquid in- 
teraction of Eq. ((83| in the form 



H 



K 



FL 



6 L 



(X0 + X')^77(X0 + X') 



(Al) 



where xo = CnR + e ^""CnL, with CnR/L the operators 
that annihilates electrons in the states with energy Cn^^ ■ 
The other modes are contained in x' = J2n'^ni^n'R + 
e~''°'Cn'L), with the sum restricted to states near the 
Fermi surface, i.e. |e„/ | ^ D, where ^ is the 
cutoff. We expand the interaction into four terms 



SHpL = SH 



FL 



SH 



(2) 
FL 



SH 



(3) 
FL 



SH 



(4) 
FL- 



(A2) 




4- 



M levels 



M levels 



FIG. 8: Occupation of the energy levels of the embedded 
quantum dot in the strong coupling limit for A*' — 4{p + n) +3 



and Of ~ 0. The pair of states with energy 
occupied by a single electron. 



is partially 



The first term, 



H 



(1) 

FL 



K 



6 L \^0 2^° 



(A3) 



only involves xo a-nd couples quasi-degenerate states with 
the same total number which differ by the distribution 
of electrons in the partially filled level. Before we look 
at the effects of this term, we argue that the other three 
terms can be neglected. The second term is 



.„(2) A<A' ,+ <? , rfO" 

^^^^ = -^ — ^ 2^ -^ 2 



X 



(A4) 



This operator is diagonal in the subspace of states with 
fixed N since it does not contain xo Sind does not act 
on the electrons occupying the partially filled level. The 
contribution to the ground state energy is 



SE 



(2) 
FL 



AttC 



K 



6 L 



:x 



(A5) 



where () denotes the expectation value in one of the states 
for fixed N. We obtain 



^41 - 



x(«!.^ + e'"<!j(VH.+ 
X (c^!„ + e*"c^I,)(c +e-*"c )\ 



K 



6 L 
3D^ 



i-X E E 1 

m— n — J\/ m^— n+1 



2T, 



K 



(A6) 



where M = DL/ {2ttvp) is the number of states above or 
below En and inside the cutoff. We choose the cutoff to 



be symmetric without loss of generality, since the effects 
of particle-hole symmetry breaking can be cast into the 
potential scattering term of Eq. ([72]) . The important 
point is that the 0(1) contribution to E{N) in Eq. (jA6P 
is cutoff dependent, but is independent of the occupation 
of the levels e^^^ (independent of £'). Therefore, it gets 
cancelled when we take the difference E{N -I- 1) — E{N) 
and has no effect on the charge steps. 
The last two terms are 



SE 



(3) 
FL 



"^FL — 



' 6 L 



Xo^X -X^^Xo 



6 



^ ^X^^Xo ■ Xli^X 



(A7) 
(A8) 



In principle, () should be regarded as a matrix in the 
subspace of quasi-degenerate states for fixed TV. Let |a)(8) 
\FS) and \b) (8) \FS) denote two states in this subspace, 
with I a) , 1 6) the states of the electrons in the partially 

filled level and |F5)^n„'<„,.=Ta. 
Fermi sea. The corresponding matrix element is 



<^'R<k |0) the filled 



SE 



(3) 
FL 



K 



6 AL 



(A9) 



Likewise, 

r(4) 



SE^;i = 



An^K 
6 4L 



iaT,i^arxMS>;(a 



At 
XoiyXo 



(AlO) 

Note that the factor of M is the same in Eqs. (jA9p 
and (jAlOp because the cutoff is particle-hole symmetric. 
Combining (|X9)) and (jAlOp . we get 



^41 + <i 



Att£_k DL 



2 

AL 



4L 2'KVF 
DS"^. 



(«|{ 



Xoa^'Xo^} &) 



(All) 



As a result, this contribution is also diagonal in the sub- 
space of fixed N and is independent of I' . Therefore, to 
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first order in ^/^/L, the splitting of tfie charge steps is 
determined by Hpj^ only. 

We now turn to the calculation of SEpj^ using degen- 
erate perturbation theory. The trivial case is the ground 
state energy for TV mod 4 = 2 (£ = 4n-|-l, n integer). 
In this case, which includes the ground state for jj, = fj,Q 
at a ~ 0, the ground state is non-degenerate: two elec- 
trons are bound in the singlet and the remaining 4(p-|-n) 
electrons fill up the single-particles energy levels of the 
ideal ring. This situation corresponds to £' = 0, and the 
contribution from 5Epj^ vanishes: 



A{p + n) + 2] = 0. 



(A12) 



Now consider N mod 4 = 3 (f = 1). For (k/L = 
0, there is one extra electron on top of the filled Fermi 
sea (as illustrated in Fig. [5]). Since the Fermi liquid 
interaction commutes with the total spin, we can work 
in the subspace where this extra electron has spin up. 
The two quasi-degenerate states the electron can occupy 
are 



W^4!h\FS) , \2)=cll\FS), 



Tt 



(A13) 



where \FS) = \GS{N = 4(p + n) + 2)) is the filled Fermi 
sea. The unperturbed energies (for ^k/L = 0) are 

En = (l\Ho\l)=E[A{p + n) + 2] + e^, (A14) 
El = {2\Ho\2)=E[4{p + n) + 2] + e^^. (A15) 



It is convenient to rewrite Hpj^ in the form 



H 



(1) 

FL 



.^^Jfc^ +e'"c^ )-(c +e-'"c ) 
6 L 

|(T |(T |(T |(T 



\{SnR + Snhf 



+2{snR + s„l) • ( e"'"cJ,H Jc„i + h.c. 



e-''"4^|c„^-cL|c„^ + /i.c.)^, (A16) 



where 



.R/L 



^nR/L 
r(l) 



(a/2)c 



nR/L 



.We calculate the ma- 



trix elements of Hpl in the subspace of states |1) , |2). 
The associated matrix for the Hamiltonian including the 
Fermi liquid interaction is 



Hn + H 



FL 



E, 



4L ^ 



4L 



Er 



4L 



(A17) 



Diagonalizing this matrix, we find that the ground state 



energy for ^' = 1 is 

E[N = A{p + n) + 3] 
= E[N = 4(p + n) + 2] + ^" 



4L 




4L 



K 



E[N = 4(p + n) + 2] + A ( 2r7. - 



A^e 



K 



4L 



-A 




(A18) 



For N mod 4 = (£' = 2), we have to distribute 



two extra electrons in the levels with energies £n^^ ■ 
£,k/L = 0, the subspace with total sfj = sj;^ + sj;^ 
is spanned by four states 



For 

= 



\RR) 
\LL) 

\RL,s) 
\RL,t) 



„it „Tt 



\FS), 
\FS), 



1 f 

1 ^ 
71 ' 



^nR^nL 
^nR'^nL 



^riR'-nL 
'^riR'-nL 



\FS), 
\FS). 



Of the above states, the first three are singlets (total 
s„ = 0) and the last one is a triplet (s„ = 1). The 
corresponding unperturbed energies are 



Err = E[N = A{p + n) + 2] + 2e^ 
Ell = E[N = A{p + n) + 2\ + 2e. 
Erl,s = ERL,t = E[N = iip + n) + 2] 



R 
n I 
L 
n ' 



There are still two other states with = ±1 that are 
quasi-degenerate with the above states. For / L ^ 0, 
the Fermi liquid interaction can only mix the states in the 
= subspace. We note that \RL, t) does not couple to 
the other three states because Hp'l commutes with the 
total spin. Moreover, it is apparent that \RL, s) can only 
couple to \RR) or \LL) via the term of Eq. (jA16p that 
transfers one electron between right and left channels, 



{SnR + Sul) 



^nR' 



h.c. 



(A19) 



However, this term involves the total spin operator s„ = 
SnR. + SnLi which annihilates singlet states. Therefore, 
\RL,s) does not couple to \RR) or \LL) either. As a re- 
sult, we just have to diagonalize the 2x2 matrix spanned 
hy {\RR) ,\LL)} 



Ho + H^pI ) 



E 



RR - -n- 

Att^k „ia 
iL ^ 



Ell-^ 



(A20) 
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The smallest eigenvalue of (jA20p yields 

E[N = 4(p + 7i) + 4] 
^ E[N ^A{p + n) + 2]+e^ + . 



K 



AL 



= E[N = 4(p + n) + 2] + 2A ( 2n + - 



Att^k 




4L 




(A21) 



The corrections to Efj^^s and iJ/jx,* &rst order in S,k I L 
are given by 



RL,s 
RL,t 



SH 
SH 



(1) 

FL 

(1) 
FL 



RL,s 



0, 



RL,t) = - 



Att£_k 
2L ' 



(A22) 
(A23) 



The ground state energy given in Eq. (jA21|) is lower than 
Erl,s, EpiL^t, except at a = 0, where it is equal to Erl^i- 
In fact, for a = 0, we have = e^. It is convenient 
to use the basis of even and odd channels, in which the 
Fermi liquid interaction of Eq. (|A16p assumes the form 



H 



(1) 

FL 



2A ^ 
3 L 



(A24) 



where s„ 



>l)(o'/4)(c„j^ + c„i) is the spin 
operator for electrons in the even channel. The ground 
state for ^' = 2 is obtained by adding one electron to the 
even channel and one electron to the odd channel. In this 
case. Sue = 1/2 and 6E = — A7r^if/(2L). For this given 
energy there are two degenerate states with = 0: 



|eo) = c}lc]^o\FS) 
\oe) - citcTt \FS) 



, Jt „Tt _„Tt 

-nR ~r ^nL ^rafl ^riL 

%/2 ^/2 

„it _„it „Tt , „Tt 

-^nR ^nL ^nR ~^ ^ 
V2 ^/2 



\FS), 
\FS). 

(A25) 



These can be recognized as linear combinations of \RL, t) 
and {\RR) + \LL))/ y/2. The latter is the eigenstate of 
the matrix in Eq. (|A20|1 with eigenvalue given by Eq. 
(jA21[) for a = 0. One can also verify that the energy of 
the states with sfj = ±1 is lowered by the same amount 
as Erl^i- For a = 0, this corresponds to putting one 
electron in the even channel and the other in odd channel, 
both with spin up or both with spin down. 

Finally, the calculation of the ground state energy for 
N mod 4 = 1 (£' = 3) is analogous to the one for N mod 



= 1). We find 

E[N = 4(p + n) + 5] 
E[N = 4(p + n) + 2] + 



3(e^ + e^) 



4L 

E[N = 4(p + n) + 2] + 3A ( 2n + I 



A< 



K 



4L 




K 



AL 



(A26) 



Now consider a « tt. In this case we have to consider 
that the level crossing involves the unperturbed levels 
e^_i(a « tt) « e^(a « tt) « 2A(n + 1/4). The trivial 
case of filled shells, in which the correction to the ground 
state energy due to the Fermi liquid interaction vanishes, 
occurs for N mod 4 = 0. In analogy with the calculation 
for a ~ 0, we find that for a « tt 



E[N = 4(39 + n) + 1] 
E[N = 4(p + 11)] + A ( 271 



A-Ki 



K 



AL 



-A 




, (A27) 



E[N = 4{p + n) + 2] 
E[N = 4(39 + n)] + 2A { 2n - 



Att^ 
4L 



-A^ 



'2{a - n)' 


1 

+ 


ft) 


TT 







,(A28) 



E{N ^ 4(p + n) + 3] 
= E\N = 4(39 + n)] + 3A ( 2n + - 



Att^ 



K 



4L 




AL 



(A29) 



The resuhs in Eqs. (Pl2|) . (^18)) . (|MT|) and (|X26)) 
apply directly to the case a ~ rmr, m even, if a is re- 
placed by a - mTT < 1. The results in Eqs (|A27IIA29|) 
apply to the case a ~ m7r, m odd, if we replace a — tt 
by a — miT. The expressions for the charge steps in Eqs. 
([5^- 157)) are obtained by taking the difference ^^+1/2 = 
E{No + £+l)-E{No + i) = EiAp + e + 2)-E{Ap + e+l). 

Note also that for ^k/L ^ the expressions for E{N) 
and the charge steps reduce to the exact ones for the 
strong coupling fixed point for all values of a. Defining 
i?^°^(a) as the ground state energy for (,k/L = 0, the 
ground state energy for ^k/L <^ 1 and a ~ mi: can be 
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written as These are the expressions used in Eqs. ((90)) and (|92p to 

calculate the persistent current. 



E{NoM) = £;^"'(iVodd) + 



4i 



[o-rHV^/^y, (A30) 



£^(7V even) even) + 2A' 



4L 



[a — rmr] \ Att^jc 



TT 4L 



2[a — m7r]\^ f Tr£,K^ ^ 



+(^) -(ASl) 
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